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ABSTRACT 


This research was concerned with the temperature dependence of thermal 
neutron flux. In view of the excellent neutron detection capability 
exhibited by the Lithium Iodide Europium enriched crystal this phosphor 
was utilized for the investigation. A miniature crystal (4mm x 4mn) 
and a flexible light pipe comination were used to determine the neutron 
spatial distribution as a function of temperature using a Pu-Be source © 
in a light-water medium. | 

A separate experimental investigation, with Indium foils and activa- 
tion analysis, was used to determine the temperature dependence of thermal 
neutron flux under the same experimental conditions as the crystal. 

A theoretical analysis was also conducted using Two-Group theory to 
obtain an expression for thermal flux with spatial and temperature 
dependent parameters. 

The results of the research indicates that the thermal neutron flux 
is temperature dependent in the region of temperature investigated 
(77° - 168° F). The experimental results verify the increasing trend 
of flux magni tude wi th increasing temperature, as predicted by the Two- 


Group solution. 








TABLE OF CONTENTS 


GMT RODUCTION. <<. anes. 1s occocematun 2. of sllcwbineus aetee Reo 


ROW, << 5 4 cs 0 0 ee 06 6 6 + OS 5 6 Semen an micemene 
Pee bHO-GROUP DU RRUSTON THEORY as so « « « + mouceunst (ool ueaeeumne 


i 
/ 


cee TEMPERATURE JDEBBNOERNCE sOG SEU TRO MAGE: sm.) ovis, « It 


C. TEMPERATURE DEPENDENCE OF SCATTERING CROSS SECTION. . 


D. TEMPERATURE DEPENDENCE OR USHeneCENGI . see 


E. TEMPERATURE DEPENDENCE OF DIFFUSION COEFFICIENT . .. . 
F, SUMMARY OF THEORETICAL ANALYSIS .... ee ee ee 


Pee DESCRIPTION OF EAPERIMENT USING DNDIUM FOILS .. . 
Pee DETEREENATION Ur SATURATED ACTIVITIES . . .. 3... 


POPMERTONAU LON GF RELAAAMIONMEERGIH . 9. 2 6 6 0 se ee 


D. DETERMINATION OF WIGRATION AREA, AGE AND DIFFUSION 
LENGTH e e e e e 6 @ e e 6 e e e e e e e e e e e e e 


Ieee eR nMEnTAL DATA WITH CRYSTAL . 2 1. ws 0 ew 8 we tw 
Rem UlocUoolON OF CRYSTAL DETECTOR . . . « «+ ce eee 


be DETERMINATION OF CRYSTAL e e@¢ ee «© © ee © @ ee @ @© @© @ @ 
Pie e tee Ole Ome APE MENG «ues. «6 4 © “66: 6 


m™m ae 7) 
® e 


» SOFT GAMSA INFLUENCE AND MULTICHANNEL ANALYZER 


PCO HeieO ram NW a ceeaie MUG? ou ie) ce! ee <s. ool esl tess 

V. COMPARISON OF THEORETICAL AND EXPERIMENTAL RESULTS .... . 
oe (neGreRiCAb VERSUS INDIUM FOIL DATA . . . ss 4 « ue 

Eee eer cn VEROUS CRYSTAL DATA %. <6 <:.s «te steer meueen 

PMO GUSMO Ges ee sk kw ee ose te ees te 


APPENDIX A = COMPUTER PROGRAM 1 FOR EVALUATION AND COMPARISON 
MME VENUES 9 es ee 8 ek tl we : 


DISCUSSIOW OF RESULTS OF CRYSTAL EXPERIMENTAL DATA. . 


13 
17 


Pale 
eie2 


26 


ee 


33] 
34 


oo 
ou 
62 
62 
63 
68 
69 


7/ 


. 84 


90 


2) 








APPENDIX B = COMPUTER PROGRAM 2 FOR THEORETICAL PLOT OF 
DIFFUSION LENGTH, AGE DIFFUSION COEFFICIENT AND 


Se IERLWG CROSS SECTION ..... 


APPENDIX C = COMPUTER PROGRAM 3 FOR THEORETICAL FLUX PLOT... . 


® @ ® a ® ® ® e ® « ® @ 


APPENDIX D = COMPUTER PROGRAM 4 FOR EVALUATION OF SATURATED 


me VITY & e ® @ ® @ @ a ® 6 ® @ e 


APPENDIX E = COMPUTER PROGRAM 5 FOR DETERMINATION OF MIGRATION 
AREA, AGE AND DIFFUSION LENGTH . 


APPERDIX F = COMPUTER PROGRAM 6 FOR FLUX PLOT FROM INDIUM FOIL 


® e ® ® @ @ ° ° ® ® « e @ e 


APPENDIX G = COMPUTER PROGRAM 7 FOR FLUX PLOT FROM CRYSTAL 


BeeePOR DATA . . « 6 6 0 3 


BMeeOGRAPHY . 2. . 2 5s 0 ee 


POITTAL DISTRIBUTION LIST . . 
FORI UD 1473 o @¢ © e@© © e€ @© @e @ 


® © e ® ® ® e ® e ® ® 


® ® e ® e ® a e e ® a ® 





ACKNOWLEDGEMENT 


A deep sense of appreciation is expressed for the continued 
guidance and inspiration given to this researcher during the period 
of investigation by my thesis advisor Profassor Dong Nguyen. Gratitude 
is also expressed to my wife Marion who provided moral support during 
periods of frustration and who assisted in the preparation of the final 


draft. 








I, INTRODUCTION 


The use of a crystal for neutron detection has normally been limited 
to large crystals which are directly connected to a photo-multiplier tube. 
Optimum performance with a crystal is reportedly obtained when the crystal 
is maintained at near liquid nitrogen temperature [1]. In that no informa- 
tion was available for usage of the crystal under the desired temperature 
range (77° to 168° F) the effect of temperature change upon the crystal's 
detection capébility had to be determined. 

The LiI(Eu) crystal as per Price [2] and Iturray [3] was first proposed 
aS a neutron detector during the nid 1950's. Because of the high 
absorption cross section, the crystal has a high detection efficiency 
and should be suitable under different experimental conditions in pulsed 
neutron research dealing with small size media and fast decaying neutron 
fields. The Boron Trifluoride filled gaseous detector can also be manu- 
factured in small sizes; however, the efficiency remains much less than 
the LiI(Eu) detector. In addition the BF. requires a minimum length of 
Signal lead between the detector and the preamplifier; whereas , the 
light pipe and crystal combination provide considerable flexibility. 

The indium foil activation analysis method is well established, but 
extremely time consuming in the ultimate determination of the thermal 
neutron flux. 

The purpose of the research was to investigate the temperature 
dependence of thermal neutron flux using the crystal detector. The 
results of the crystal data are to be compared with a theoretical 
analysis and additional experimental data obtained by the indium foil 


activation method. 








II, THEORY 


A, TWO GROUP DIFFUSION THEORY 


The object of this theoretical analysis was to develop an expression 


for neutron flux with temperature dependent parameters for comparison 
with data obtained from experimental phases of the investigation. The 
analysis considers a point source in an infinite medium of water. 
Neutrons (are born) at the source with high energies. The average 
neutron enercy is 4.3 Kiev for Pu-Be source. As a result of collisions 
with the water molecules the neutrons are either absorbed or moderated 
with an enercy loss. The steady-state neutron flux was investigated 
with energy end position as variables. A two group metnod of analysis 
utilizing the steady-state diffusion equation was used to develop an 
expression for the neutron flux. Group 1 neutrons (Fast Neutrons) are 
categorically those neutrons with energies greater than | ev and Group 
neutrons (Thermal Neutrons) are those neutrons with energies less than 
lev. The diffusion equation for Group 1 can be expressed as follows 
where D is the diffusion coefficient, De is the total scattering cross 


section, and ¢ is the flux with energy and position dependence. 


Dy(E) V6, (rsE) = ECE) $y (rs) 


{ E(E' > E)d,(r,E')dE* = 0 
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Equation (1) can be simplified by assuming that the flux is separable 


in space and energy and replacing o(r,E) by F(r)X(E). Integrating al} 
the terms in equation (1) with respect to energy over the limits of 


the Group 1] neutrons results in the following: 
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Substitution of the aforementioned relationship into equation (2) results 
in the following differential equation: 
2 Tene eye 
Div Fi(r) - La Fy(r) = 0 


or 
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Ty being the age of fast neutrons to thermal energy. We now return to 
the more conventional symbol for the spatial neutron flux and write 


the above equation as: 
v64(r) - “oy Cr) =a (3) 


The general solution to equation (3) in spherical coordinates is as 


follows: 
“Kir kk 


(ry) = + (4) 


The constants can be evaluated with the following boundary conditions: 


7 $4(7) must be finite as r tends to ; therefore, C, = Oe 


2. the limfar“d, (1) = §, where J = -D grad oy and S = source 
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In that the thermal flux is dependent upon the fast flux, a solution 


for the thermal flux $5 (1) was assumed and then verified. 
2 2 > 
V dol) = Ko ¢o(r) + D, $4 (1) = Q (6) 


¢o(t) = S3o4(r) + G2Zo(r) (7) 





Substitution of equation (7) into equation (6) and collecting similar 


terms results in the following expression: 


| b 
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By comparison of the coefficients of the fest flux of equations (3) and 
(8), the coupling coefficient (S3) must satisfy the following relation 
Ship: 


d D,/2 
2 enc eles 
ot ia Fe Z 


2 Trek 


fecnat Ko 


2. 2p? a I, 
L 
2 


on 1, 
v Zo(r) - ky 


Z,(¥) = Q 


The solution of the latter equation in spherical coordinates is 
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The boundary conditions for thermal flux for evaluation of the constants 
are as follows: 


il. o5(r) must be finite except at origin; therefore, Co = Oe 


a do(r) for a point source is zero at the point source 


; Z 
1,.6., lim 4cr Joli); where Jo(r) = -D, grad $5(1) 


ye 0 





The final solution is: 


Z er/L  -r// T 
SL fe - ag, 7] 


(Lars) (9) 
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The physical parameters +, L, and D are now investigated for their 


temperature depenaence. 


Be TEMPERATURE DEPENDENCE OF NEUTRON AGE 

Neutron age does not represent an clapsed tine, but is proportional 
to the mean squared distance tnat a neutron travels between formation 
at birth with high energy (at the source) until thermal energy is obtained 
(approximately 0.025 ev) as a result of being moderated by the medium 
Surrounding tne source. A small value of age indicates that the neutrons 
have suffered little slowing down and have not diffused far from the 
Source. The magnitude for age in a water moderator for Pu-Be source 
is given as 52.8 em? [5] at ambient temperature. The temperature 


dependence of age is developed using the following basic relationships: 
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Taking the natural log and differentiating with respect to temperature, 


it is obtained d( ) 
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where N is the numper of moderator nuclei per unit volume and Oo. and 
o,, are the microscopic Scattering and transport cross sections 


respectively. Assuming that O. and o,., are only slightly temperature 


tr 
dependent in the range of interest the following approximation is made: 


nv (10) 


Since the numer of Nuclei (HN) is inversely proportional to the volume, 


it follows that 


or , letting a, be the temperature coefficient of volume expansion, 
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Over the region of temperature range investigated (25° C to about 75° C) 
the temperature coefficient of expansion (a). is nearly linear [4] and 


1S approximated by the following expression: 


z ets met " 
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Substituting the expression for the temperature coefficient into 


equation (10) and integrating yields the following: 


<(T) = C EXP [1.8 x 107° T + 1.04 x 107°T*] 
Using the experimental value of age to indium resonance for a Pu-Be 
source (at T = 25 C) the constant was evaluated and the final expression 


for age with temperature dependence was determined: 


We 





7(25° C) = 52.8 + 2.5 cm* 
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c(T) = 52.2 EXP [1.8 x 107'T + 1.04 x 107°T*] (12) 


C. TEMPERATURE DEPENDENCE OF SCATTERING CROSS SECTION 

Prior to the investigation of the temperature dependence, the energy 
dependence of the macroscopic scattering cross section ie) Was 
considered. The microscopic scattering cross section of water decreases 
witn increasing energy. The elastic scattering cross section is not 
usually measured directly because of experimental difficulties and is 
usually obtained by subtracting the cross section of all other processes 
(capture, inelastic, etc.) from the total cross section. The total 
cross section of water is essentially all scattering. The following 
empirical relation has been proposed [4] to express the scattering cross 


section as a function of energy in a water moderator. 
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An analysis with a Taylor series expansion along with an iteration 
technique using the conputer was conducted to check the non-integer m. 


Experimental data for microscopic scattering cross section! were fitted 


| Approximate values for microscopic scattering cross section as a function of 
eneray were obtained from a curve contained in BNL-325 [7] by taking the 
difference of the total cross section (o-) and the absorption cross section 


iS 


by assuming an expression of similar form to the given empirical rela- 
tion. A relationship was obtained by this procedure as shown in Appendix A. 
In order to attain greater accuracy the energy dependent scattering 

cross section data from the computer program Gasket-Flange II [6] was 

used and similar results were obtained. Two triple curve plots (Figs. 1 
and 2) were plotted by the computer. One curve was based upon the 

FLANGE II data which represents the microscopic scattering cross section 
as a function of energy based upon the best available information. The . 
second curve is dependent upon Sees relation (i = 0.281) obtained | 
by fitting measured scattering data (0 to .5 ev range) for water. 

Figure 1 shows that the value of m = 0.281 yields results very close 

to the Gasket-Flange II data, The validity of the empirical non-integer 
(mm) of Ref. 1 was not confirmed. Further investigation was conducted 

by using extended energy ranges (Fig. 2) up to 1.5 ev and a revised m 
(,255); however, similar results were obtained. All functions that are 
dependent upon the non-integer (m) were evaluated, using both values of m. 
The temperature dependence of the scattering cross section was analyzed 


in a similar manner to that used for age: 
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where O. was established in Section B to be 


a = 9.0 x 107° + 1.04 x 107° T 
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Substituting the temperature coefficient into equation (13) and integrating, 


E(T) = C EXP [-9.0 x 10 T - 5.2 x 10° 1) 


The constant is evaluated at thermal energy with the following values: 


] 


r= 3.443 cm Eee) 00025 even =a2028G) 


So isG eee Or 10m eno Ee ee 
S e 4A e \ ¢ (14) 


The macroscopic cross section is illustrated in Fig. 3 as a function of 


temperature for thermal energy. 


D. TEMPERATURE DEPENDENCE OF DIFFUSION LENGTH 

The diffusion length is one-sixth of the mean square aistance that 
a neutron travels from the point at which it becomes thermal to the point 
at wnich it is captured. 

An empirical relationship based upon prior experimental observations [5] 
with a Pu-Be source in a container of similar geometry was used to express 


tne temperature dependence of the diffusion length as follows: 


L(T) = 2.64 + 0.0061 T cm | (15) 


where T 18 in decrees centigrade. 


E. TEMPERATURE DEPENDENCE OF DIFFUSION COEFFICIENT 
The temperature dependence of the diffusion coefficient for the 


thermal flux was studied by using the following empirical relation [4]: 


D(E) = D(E,)(E)" 
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By definition 
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Assuming a Maxwellian energy spectrum: 
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For material having essentially constant scattering cross section, 

m= 0 and the average thermal diffusion coefficient is temperature 
independent; hovever, for a water medium the non-integer (ra) is not 

equal to zero and the diffusion coefficient is temperature dependent. 

As outlined in Section C, the empirical value for m was 0.47 and the value 
determined by fitting data from FLANGE II was determined to be 0.281. 

By normalizing the diffusion coefficient to the value at the reference 
temperature To = 25° C, a temperature dependent expression is obtained 


as follows: 


D(To) = .16 cm [4] 
_ . 16 
D(Eo) ~ F(m + 2) 
D(Eo) = 0.1229 cm m= .47 


D(Eo) = 0.1387 cm m= .231 
A graphical presentation of the diffusion length as a function of temperature 


is illustrated in Fiy. 4. 6 
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The temperature dependence of the diffusion coefficient is illustrated 


aeeiG. 5. 


Fo SUMMARY OF THEORETICAL ANALYSIS 

The two-group diffusion theory was used to obtain an expression for 
the thermal flux. The parameters of this solution were investigated 
for their temperature dependence. Several computer programs (Appendices 
A, B, and C) were used for compilation and tabulation of data along 
with providing graphical presentations of all parameters as a function 
of temperature. 

The graphical presentation of the computer plotted thermal flux 
as a function of distance from the source and temperature illustrates 


only a slight increase in flux with an increase tn temperature (Figs. 6 


and 7). 
An additional illustration relating the theoretical temperature 


dependence of related parameters; namely, age, diffusion length and 


2 2 


migration area (nie), where [7 = ty ee cme. is contained witnin Fig. 8. 
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IIL. EXPERIMENTAL DATA WITH INDIUM FOILS 


A. DESCRIPTION OF EXPERIMENT USING INDIUM FOILS 

A VISIFLUX neutron howitzer was utilized to conduct the flux 
measurement us ing the Plutonium-Berylium (Pu-Be) source and Indium 
foils. The VISIFLUX howitzer is a plexiglass cylindrical tank having 
a2 foot diameter and a 2 174 foot height and is specifically designed 
for nuclear experimentation. A water level of 45 gallons was maintained 
during the periods of foil irradiations. [Electric ‘craton type heaters 
along with a control unit were used to obtain and stabilize the desired 
operating tenperature. An agitation type mixer was used to insure uniform 
tank temperature. The VISIFLUX howitzer and associated equipment are 
miustrated in Fig. 9. 

The Indium foils were 1.0 inch in diameter and 0.010 inch thick. 

A total of 20 foils were used with 10 foils enclosed within individual 
cadmium covers of 0.020 inches thickness. Plastic foil holders were 
used to place the bare foils on one side of the source and the cadmium- 
covered foils on the opposite side of the source in a horizontal plane 
perpendicular to the center of the source. The foil Teen were used 
as criteria for soa Hon ae with the lightest foil being placed nearest 
the source. 

Both sets of foils were irradiated at the same time at a stabilized, 
predetermined temperature for a minimum period of 3 hours. The activity 
was determined by using an end-window Geiger-Mueller tube in a vertical 
lead shielded housing with an electronic timer-scaler counting device. 
Specific equipment used for the Indium foil experiment is listed in 


Table 1 and illustrated in Fig. 10. 
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TABLE 1 


EQUIPMENT 
ITE‘ MANUFACTURER MODEL 
Timer-Scaler Hewlett-Packard 5201L 
Preamplifier Hewlett-Packard 55544 
GM tube Amperex 200HB 
Vertical Lead Shield Designers for Industry AL=12A 


High Voltage Power Supply Power Designs Pacific, Inc.  2k-10 


The same procedure was adhered to for each of the measuremants at 
different temoeratures. Background readings were obtained randomly 
before and during periods of irradiation. Background counts were found 
to vary considerably during the day time while in general, background counts 
were lower during the evening hours. The procedure was to count the 
bare foils first, starting with the foil located closest to the source. 
Each foil was allowed to decay for three minutes to insure the desired 
count rate from the 54 minute activity was not influenced by the 13 
second half life of tn. To attain maximum consistency in the counting 
procedure all foils were counted for five minutes with the surface of 

the foil closest to the source being placed facing the G-M tube. 

Prior to conducting the experiment a preliminary investigation was 
conducted on the counting equipment to determine the G-l tube plateau 

for proper operating voltage, and dead time correction factor for the 
counting equipment and also a Statistical analysis to establish the 


presence of other than statistical variations. 
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The operating voltage for the G-M tube was determined by plotting 
observed count-rate versus increasing high-voltage [2]. The G-IM 
plateau voltage was determined to be 1200 volts. 

A dead-time correction factor is required in that the period 
immediately fcllowing a G-l tube discharge the electric field is 
reduced because of presence of positive ions. Any particle entering 
the tube during the pericd of reduced electric field will not trigger 
a discharge and subsequently the observed counting rate will be less 
mien the actual count rate. A split-source of radium was used to 
determine the dead-time correction factor [2] with a period of 38 micro- 
seconds. To verify that the counting equipment was operating normally, 
the following technique was applied [2]. 


The standard deviation is defined by the expressions: 
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ry and ro are observed count rates using a co’ (0.97 
' micro-curie) source. 
ty and t, -are length of time of observation . 

The method consists of calculating the standard error of the 
difference between two counting rates which are obtained from the same 
source under identical circumstances. The standard error is then 
compared with the observed difference between the two counting rates. 

A probability table (Table 2) is used to determine if the deviations 
are due to causes other than the statistical nature of the disintegration 


process. 
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TABLE 2 
PROBABILITY P OF OBTAINING /Y, “Yo > ko 


K eu Wee 1.4 1.6 240 oa 3.0 


P OF 159 0 ns 0.03] 0.055 R023 0.006 0.001 


The results of the analysis indicate that the counting equipment was 
not subject to receiving erratic counts and implies the error incurred 
will only be of a statistical nature. This technique was applied 
several times during the course of experimentation with similar results 


obtained, 


B. DETERMINATION OF SATURATED ACTIVITIES AT DIFFERENT TEMPERATURES 
The purpose of Indium foil activation was to determine the neutron 
distribution as a function of both temperature and position. The 


M6, 


production rate of the Indium isotope (In _is proportional to the 


neutron flux density and the activation cross section, 7.¢e., 


P= zo oV 

where P = production rate 
aes macroscopic activation cross section of Indium 
@ = total neutron flux density 


V = volume of Indium foi] 


In that the isotope (in! |) 


has a 54 minute half-life and is subject 
to decay wnile being irradiated the following expression represents the 
production process 


dN _ . 
at = 2a¢¥ = AN (18) 


ai 








where 


dN 
di 


i 


production rate 


AN = rate of decay 
4 = decay constant (1.28 x 1072 min?!) 
Integration of equation (18) yields the following: 


: =r¢ 
= z, oVO oe) 


] r 
where 
Ny = the no. of atoms of Int !6 existing at time t, 
AN, = the absolute activity at time t 
Sue that 


Ay = kxily = ku, V(1-e7"1) 
where k is a proportionality constant between the absolute activity (4) 
and the counting activity (A,). 
For a large irradiation time ty the saturated activity would be 
obtained, i.e., AS = KE OV Such that the activity at time ty can be 


represented by A, = As(1 - acee 


In that the saturated activity is 
the quantity being sought and Ay is the actual extent of activation we 


express the saturated activity in terms of the actual activity: 


As = AV/(1 = 7") 
As mentioned previously, a waiting period of about three minutes was 
applied because of a Beta emission with a 14 second half-life. At the 
time the counting begins (t,) the activity of the foil has been reduced 
and is represented by 


: — arat,, “A(t, - t,) 
A, = AA] enn je 2 1 


Sc 








such that 
Ao 
ALF 
S (] sf an Abr) (7 Alts - ty)) 


Noting that A, is the activity at the time ty and that the actual count 
rate as observed by the detector is for the period (t, ‘ to)s the total 


count (C) can be calculated as follows: 


Coane 


3 2 
C | A(t) dt 


0 
where ty - te = counting time 


C = observed activation count 


het) = activity 


C. ~ t, 
pr CO = AC - a7 Aty yen Al ty i ty) [ors 
0 
such that 
WE 
ar | ae (19) 
eee KC? Taare 3 1) 
where 


t = {rradiation tine 


to ty decay tiie 


ce see 


1 decay time + counting time 


The same expression was used for the cadmium-covered foils. A correction 
factor for the absorption of fast neutrons by the cadmium cover based 
on foil thickness was obtained from the interpolation of a graph by Tittle 
as presented by Price [2]. 
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The observed count rates (corrected for background) as a fundtion 
of distance were used as input data for a computer program listed in 
Appendix D. ‘This computer program provided the following data for both 
the bare and cadmium covered foils as a function of distance: 


1) saturated activity 


NO 


3 


) 

) saturated activity times distance squared 

) saturated activity tines distande to the fourth power 
) 


4) the standard deviation of the saturated activity which is 
designated as sigma (oc) | 

Four graphical presentations were obtained: a semi-logarithmic plot 
of both the bare and cadmium-covered activities times distance squared 


and linear plots of the same curves for ultimate area integration. 


C. DETERMINATION OF RELAXATION LENGTH WITH TEMPERATURE DEPENDENCE 

In contrast to charged particles which have a defined range, neutrons 
and gamma rays are stopped approximately in an exponential manner. The 
relaxation length is the thickness of a material that causes a drop in 
intensity by a factor of e owing to slowing down or absorption [8]. 
The relaxation length is energy dependent and will have a specific 
length for a specific known neutron energy; namely, the relaxation 
length of thermal neutrons (E = 0.025 eV) is equal to the thermal 
diffusion length and for Pu-Be neutrons (E = 4MeV) the relaxation length 
is of the order of 10 ci. 

In view of the nearly linear temperature-dependent behavior of the 
theoretical relationship established for the macroscopic scattering 


cross section (Ref. II-C) and the inverse relationship of the relaxation 
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length, b = 1/zs(E), it was assumed that the relaxation length would be 
essentially linear with temperature over the range of temperatures 
inves tigated, 

At large distances from the source the activity curve is exponential 
in nature. Tne spatial dependence of the slowing down density (q) 
resulting fron a point source (at large distance from the source) is 


->sV 
given by q = = » In that the thermal flux is proportional to the 
r 


Slowing down density the following relationship was obtained. 


Pee caer 


or enis(r)r° = -Dsr 


and ends(r)r? x =P/b 


where b is the relaxation length, and As(r) is the activation activity. 
The problem that remained to be resolved was the determination of the 


value r beyond which equation (2) became valid. 


D. DETERMINATION OF MIGRATION AREA, AGE AND DIFFUSION LENGTH 

The Fermi Age (t(E,,E)) is defined as being equal to one-sixth of 
the average crow flight distance from the point source where a neutron 
is born (E.) to the point where the neutron acquires the energy E [4]. 
A similar relation holds true for the Migration Area except the distance 
extends to the point of capture. The difference between iligration 
Area and Age is the diffusion length squared. Relationships between 
area, lengths and mean square distances are listed as follows: 


t= 1/6 re c = age Yc" = slowing-down length 





We = 1/6 ¥ Mo = Migration area, M = Nigration length 
Co ~2 aes . : ' 

L~ = 1/6 ren L = Diffusion length 

Wee t+ 15 

ee 

re = 0s * Pep 


-2 ; hs 
The mean square distances ry can be expressed as a function of activity 


by the following relationship: 


oo 


, | 
[ Asy dr 
re =O _______ = gif 





Asr dr 


oa 


A similar relationship holds true for the age using the cadmium activity 


data such that 


co 


{ As(cd)r“dr 
re ean ne = = OT 
[ As(cd) r@dr 
0 
To account for the eneray difference (1.45 to 0.025 eV) and neutron 
travel between the Indium resonance and thermal energy a value of 1.8 cm? 
is added to the computed value of age as recomnended by Valente [9] such 


that 


L = (t + 1.8) 


The computation of uf, t, and L were conducted by the computer (Appendix £). 
To evaluate the moments of activities (Asr“, Asr’, As(cd)r“, and As(cd)r*) 


the area under the curves was determined by a Simpson's approximation. 


36 





Two segments of the curve were considered in eValuating the area; namely, 
the experimental data witn a smooth curve crawn through the data points 
and the exponential portion which is proportional to eT K/D An intercept 
and slope were determined from a computer semi-locarithmic plot of 
activity times distance Squared versus distance and used to construct 


an expression of the form ce" F/0 


to represent the exponential portion of 
the curve in tne integration process. The constant C was actually 
expressed as an exponential in view of the magnitudes incurred (Appendices D 


and E are gerinane). 


ee DISCUSSION OF EXPERIMENTAL RESULTS 

Five sets of experimental observations at different temperatures 
(77, 102, 129, 145, 168°F) were obtained and analyzed. Three graphical 
presentations for each run observed (Figs. 11 through 25) in addition to 
tabulated data (Appendix D) are presented. All the graphical plots and 
analyzed data were obtained by computer programming (Appendices D, €). 
Figures 11 through 25, in addition to the experimental data points, 
also have a smooth curve. The smooth curve is actually comprised of 
two curves, one light and one heavy. The light curve was determined by 
the relaxation length parameter (b) and was superimposed by the computer 
on the experimental data points as an aid in resolving the exponential 
curve fit. The heavier curve was drawn through the data points to optain 
magnitude versus distance such that an area integration using Simpson's 
approximation with the computer could be accomplished. The results of 
the evaluation of the area moments for ultimate determination of the 
migration area, age and diffusion lengths are listed in Table 3 and 
reflect an increasing trend with tenperature. The relaxation lengths 


are also included within the sane table. 
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An error analysis based upon the statistical error (sigma) was 
conducted for one set of experimental data (T = 102°F) and is illustrated 
in Fig. 26. Two points of the experimental data do not fall within 
the region where the smooth curves have been drawn. These points were 
assumed to reflect poor experimental results and were not considered in 
the area analysis. 

The additional errors considered were inaccuracies of foil positionirg 
(40.1 inch), foil removal/cecay timing (£5 sec), temperature variations | 
(+1°F) and human error. 

An additional analysis of the Indium foil data was conducted to 
determine the thermal neutron flux. The neutron flux was obtained from 


the saturated activities in the following manner: 


(r) = CAs(r) 

where e(r) reflects the total neutron flux 
C is a proportionality constant dependent on the Indium foils 
As is the saturated activity 


The proportionality constant was evaluated by the following expression: 


Q=4nzal{r)r‘dr = annac{ As (r) ran [n/min] 
0 O 
or ¢ = ——t__ [ofa] 
inzal As(r)r“dr 
0 
ae: anc - | 
where Bea = lece a VU cm 
Fa = macroscopic absorption cross section for water rem” * 


Q = neutron emission from source [neutrons/min J 
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ERROR ANALYSIS BASED ON: 
STANDARD. DEVIATION USING 
MIGRATION AREA CURVE (T = 102° F) 
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As(v) “dr =i2el 12x 10/ [count+cm?/min] as determined from 
computer data (Appendix E) for T = 102°F 
Such that 


8 
Hace x_10 ane Aan 


C= 9 7 
An x 2yec Oe? 1? st 





A computer program (Appendix F) was used to evaluate the total and 
Fast fluxes from their respective activities as a function of distance 
and temperature and also provide tabulated data (Table 4) and oraphical 
presentation (Figs. 27-29). Computer computations also included the 
thermal flux by taking the difference of the total and the fast fluxes, 
The graphical plots of the fluxes do not clearly show the increasing trend 
of flux magnitude with temperature because the temperature increase 
is small (range of temperature investigated is small). In order to draw 
conclusions about the measured foil data, the tabulated data as presented 
in Table 4 must be analyzed. A review of this data reflects the general 
trend of increasing flux magnitude as a function of temperature. The 
experimental values of relaxation length as listed in Table 3 were 
Investigated for temperature dependence. Figure 30 illustrates the 
linear relationship with temperature dependence which was found to be 
prevalent for relaxation length. Two curves are contained within the 
illustration and represent the relaxation lengths of the bare foils 
and the cadmium-covered foils. The slopes of both curves were found to 
be approximately the same. The difference in magnitudes was attributed 
to the smaller neutron population as seen by the cadmium=covered foil 
in that the relaxation length is inversely proportional to the scattering 


cross section of the surrounding medium. The assumption was made that 
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the relaxation length would be essentially linear in the temperature 
range investigated and this assumption is confirmed by the experimental 
results. Correlation of the theoretical versus the experimental data 


Will be discussed in the comparison of results. 
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IV. “EXPERUMENTAL DATA WITH CRYSTAL 


A, DISCUSSION OF CRYSTAL DETECTOR 

The basis of a scintillation counting system is the inherant 
characteristic of a pnosphor to convert into light some fraction of 
the energy lost by jonization of a charged particle, while the particle 
passes through the phosphor material. The light that is emitted is 
detected by a photomultiplier tube (Pii) and conmerted into an electrical 
pulse which 1S proportional to the manutd and duration of the light | 
Source. The phosphor (crystal) used for the experiment was Lithium-lodide 
with Europium enrichment (Lil(Eu)). 

The interaction of the phosphor with the neutrons is as follows [1]: 

“Lio + nt > oHet + GH + 

To attain a greater freedom of manipulation in the source of experimenta- 
tion a 5-1/2 foot section of flexible 1/4 inch light pipe was used to 
transmit the light signal from the crystal to the photomultiplier tube. 
An aluminum holder was used to provide a means of securing the light pipe 
to the photomultiplier tube. A light film of silicon crease was placed 
between the light pipe and the Pii tube to insure maximum surface 
contact. The Pil tube was then sealed by wrapping the tube with black 
rubber sheeting to insure no light leakage aside from the signal 
source received via the light pipe. The light pipe was in turn 
connected directly to the crystal by means of a hollow cylinder. 
The light pipe was secured to the aluminum cylinder by set screws. 
The crystal was placed inside the cylinder flush-up against the 


light pipe and was eld secure by a washer, spring 


G2 


and screw arrangement. Silicon grease was placed between adjoining 
surfaces of the light pipe and crystal. The arrangement describea above 
provided easy access to the crystal which was required for insertion 

of a cadmium cover when measuring fast activities. The cadmium cylinder 
and washer were fabricated to fit around the crystal. The specific 
equipment used in this experiment is summarized in Table 5 and illustrated 


mersgs. 31 and 32. 


igre, o 
CRYSTAL UEMECTOR EGUDPENT 


SO ee ae EE Ss ee ae. 


ITEM DIMENSION MANUFACTURER MODEL 
Crystal A am Harshaw Chemical Lil (Eu) 
Light Pipe 5-1/2" x 1/4" = American Optical Co. --- 
Photomultiplier Dumont 6292 
Preamplifier Hewlett-Packard 10615A 
Voltage Supply Hewlett-Packard 6516A 
Timer-Scaler: Hewlett-Packard 5201L 
Chart Recorder aa Muse ley 71008 
Pu-Be Source (3 curie) Baird-Atomic, Inc. NUMEC-M 


B. DETERMINATION OF CRYSTAL TEMPERATURE DEPENDENCE 

The Lil(Eu) crystal's performance has been shown to be optimum at 
liquid nitrogen temperature [1]. The question that had to be resolved 
Was to determine to wnat extent, if any, a variation in the temperature 


of the crystal affected the counting efficiency of the crystal. 
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To establish the nature of temperature dependence on the crystal 
the following system was fabricated. Approximately 10 feet of flexible 
plastic (TYGON) tubing (1-1/4 inch OD) was bent in the shape of a U and 
partially inserted into the neutron howitzer, A rubber stopper was 
placed on the light pipe between the crystal and the PM tube. A thermo- 
couple was inserted within the cylindrical aluminum holder containing 
the crystal such that the crystal temperature could readily be observed 
on a pen recorder. The aluminum holder containing the crystal and thermo- 
couple and attached to the light pipe was inserted into the TYGON tubing 
and placed in close proximity to the source. The stopper was inserted 
into the end of the tubing and sealed. A Tee connection mating a smaller 
section of flexible tubing (1/2 inch 0D) to the larger tubing was installed 
just below the rubber stopper sealing the light pipe. 

The 1/2 inch section of tubing was connected to a water faucet. 
The remaining end of tne | 1/4 inch section was returned to the sink 
to function as a return line. By using hot and cold water the temperature 
in the vicinity of the crystal could be varied over a moderate range 
(63° F to 131° F) corresponding to the region of temperature interest. 
Figure 33 illustrates the arrangement of the crystal and the Flexible 
tubing utilized to determine the effect of temperature on tne crystal 
counting efficiency. The VISIFLUX howitzer was filled with 45 gallons 
of water which was maintained at a constant temperature. A Pu-Be source 
was placed within the source holder after the crystal was positioned 
adjacent midway between the source and the edge of the tank. 

The aforementioned arrangement allowed a series of observations 
to be made at a constant neutron flux while the actual crystal temperature 


was varied over a noderate range. 
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The observations indicated only random minor differences which in 
view of the large count rate are attributed to statistical error only. 
The conclusion drawn from this phase of the investigation was that the 
crystal counting characteristic is independent of temperature for the 
range of temperatures investigated, Experimental data as presented by 
Murray [3] partly substantiates the aforementioned; hovever, most 


of Murray's research pertains to temperatures well below zero. 


See besCRIPTION OF THE EXPERIMENT | 

The major advantage of the phosphor detector -and light pipe .arrange- 
Ment was the capability of obtaining considerable data in a short period 
of time as compared to the foil activation method of experimentation. 

Measurement of neutron fluxes was obtained with the crystal detector 
by positioning the detector at varying distances from the source. The 
detector was shifted in a Honontat plane perpendicular to the center 
of the source. For comparison purposes the conditions at which foil 
measurements were carried out (position, water level, etc.) were duplicated 
as closely as possible. 

The detection equipment was carefully tested for proper operating 
condition in the same manner as used for the Indium foil detection 
equipment; the results indicated that only statistical variations were 
in evidence. Two separate sets of readings of count rate versus distance 
were obtained for each temperature investigated. One set of readings 
reflected the total neutron population as seen by the crystal; whereas, 
the second set of readings represented only the fast neutron population 
as seen by the cadmium covered crystal. Background readings were taken 


at random before, during and after count rate observations. The recorded 
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data was corrected for background and the thermal flux was determined by 
taking the difference between corrected total and fast observed experi- 
mental values as listed in Tables 6 and 7. Table 8 lists the thermal 
flux with distance from the source in centimeters for comparison purposes 


with thermal fluxes obtained by alternate methods. 


D. DISCUSSION OF RESULTS OF CRYSTAL EXPERIMENTAL DATA 

A computer program was written for evaluating and plotting of 
the total, fast and thermal flux magnitudes versus distance at differeni 
temperatures. fhe computer program and experimental observations are 
listed in Appendix C and shown in Figs. 34, 35 and 36. The experimental 
data indicates a general increase in nriagnitude with an increase in 
temperature. The observed aata as listed in Tables 6 and 7 reflect 
some inconsistencies in indicated counting rates; namely, some activities 
were found to be smaller in magnitude for the next higher temperature 
investigated. The aforementioned inconsistencies are attributed to a 
combination of statistical errors and positioning errors (location of 
detector from source); however, the overall trend of increased count 
rate as a function of temperature is apparent and consistent witn the 
Two-Group solution as discussed in Section II-F. | 

The experimental data is only relative in magnitude in that an 
absolute monoenergetic neutron source was not available for calibration 
purposes. Collimated neutrons when incident nornal to the crystal 
surface are absorbed at a rate proportional to the number of absorbing 
atoms present. The following expression relates the incident neutrons 
(Io) to the number of neutrons transmitted (I) 


a (18) 
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FAST FLUX AS A FUNCTION OF 
DISTANCE AND TEMPERATURE AS 
DETTE RMTAED WITH CADMIUI-COVERED 
CRYSTAL DETECTOR 
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TABLE 8 
CRYSTAL THERMAL FLUX HAGHITUDES 
AS A FUNCTION OF DISTANCE AND TEMPERATURE 


Distance Temperature (°F) 

cm 34 102 129 145 168 
6.349 975520 1023784 1025136 +=: 1020432.~=—Ss«4027832 
7.619 712392 813334 856076 876076 893432 
8,389 579838 680748 727124 699516 739672 
10.16 423680 554124. 573552 557172 595976 
11.42 379072 428048 444364 443952 476740 
12.70 246504 345052 343024 349240 375272 
13.96 236824 253008 268632 271804 289756 
15.24 172904 190492 194924 210760 221424 
16.50 147192 148500 152712 166132 174580 
ie 77 112080 112092 113700 128020 134404 
19.04 74584 86700 87372 100596 107528 
20.3) 58360 63952 69904 78796 81572 
21.58 49112 52620 52376 61336 66140 
22.85 39224 42320 45408 47820 50724 
24,12 30068 ~ 35044 33828 38460 40308 
25,39 23552 27000 24712 30784 31356 


i 


Zi = linear absorption coefficient 
X = thickness of crystal 
The crystal used during this phase of tlie investigation was a 
hermetically sealed unit in that the Lil(Eu) phosphor is extremely 
hygroscopic (moisture sensitive). The crystal thickness (X) was on 
_ the order of 4 millimeters. The linear absorption coefficient is reported 
as 16.7 cm” | [1]. The composition of the enriched Lil crystal along 


with relative properties [1] are listed in Table 9. 


TABLE 9 
Relative | Ni. oi a 
Nucleus Concentration (n/cm™ ) (barns) (cm) 
ny 4.4 8x10" 0.033 x 
he 95.62 | 1.74x10°4 950 16.5 
100.0 1..82x10°¢ Bee 0.115 
Eu 0.05 1.09x10 7 4300 0.048 
16.663 


The efficiency for the LiI(Eu) crystal can be computed by the 
following expression [1]: 


_ (Noi ti? 


as [1 - ere 


R 


ley a 16.7(.4)- — O, 
n iGo [1 C : Cai 


The calculated value for the efficiency is impressive; however, the 
size of the crystal must be considered. The number of neutrons that 
are transmitted through the crystal (do not interact) can be determined 
by evaluating equation (18) as follows: 


Teceromens > \°=") = § 00125 To 
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The number of neutrons that interact = (1 - 0.00125)Io, thus further 
enphasizing the calculaticn for high efficiency. 

The efficiency equation applies rigorously to a crystal with a very 
large radius. Since tne crystal used for this phase of the experiment 
had a very small radius (4mm) it is expected that many neutrons escape 


the crystal radially resulting in a reduced efficiency. 


E., SOFT GAMMA INFLUENCE AND WULTICHANHEL ANALYZER ANALYSIS OF DATA 

The Pu-Be source used for the experimental phase of the investigation 
was a 3 curie source with an approximate yield of 5.4x10° neutrons per 
second and a half-life of 24,400 years distributed by Baird-Atomic, Inc. 
(Model WUMEC-M). The neutrons are obtained from the Pu-Be source by 


the following reactions 


Co 


gqbu 7 lie 


Be? + 


f 2 


Energetic Gamma rays can also react with beryllium to emit neutrons; 
however, the gamna rays of the plutonium do not have adequate energy 
for the reaction to occur and are distinguished as soft gamma rays. 

In an attempt to ascertain the effectiveness of the cadmium-covered 
crystal a multichannel analyzer was used to obtain the spectrum of 

the crystal. Figure 37 illustrates the results obtained on the multi- 
channel analyzer and represents two separate input presentations; 
namely, the crystal without cadmium and the crystal with cadmium. The 
spectruin results indicate that the cadmium cover has removed virtually 


all the thermal neutrons and that the crystal is subjected to only fast 


neutrons when enclosed by the cadmium, 
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A discrimination technique wnereby the ganma influence can be 
separated from the neutron population as proposed by Furuta, Kinbara, 
and Kaieda [10] was considered for a more sophisticated approach to the 
analysis. This technique was not applied and the data presented merely 
reflects the trend of the influence of the cadmium cover with the crystal 


as observed on the spectrum analyzer. 


12 


Ve COMPARISON OF THEORETICAL AND EXPERIMENTAL RESULTS 


A. THEORETICAL VERSUS INDIUM FOIL DATA 

The theoretical and the Indium foil data show excellent correlation. 
in comparisoa of flux magnitudes. A pair of curves (T=102° F) represent- 
ing the thermal flux as determined from theory and by Indium experimental 
data are illustrated in Fig, 38. The curves were plotted from the 
computer data (Appendices C and F) as listed in Tables 4 and 10. This 
tabulated data has been rounded off for simplicity of presentation. 

The temperature trend of both the theoretical and the foil thermal 
fluxes indicate an increase in magnitude at higher temperatures. 

The experimental values of migration area, thermal age and diffusio‘ 
length have deen superimposed upon tne theoretical curves of these 
temperature dependent parameters as per Figs. 39 and 40. The experimental 
values with the exception of diffusion length were found to be sligntly 
higher than the theoretical values. The discrepancies between the 
theoretical and experimental data are mostly attributed to one assuimp- 
tion used in the derivation of the theoretical expression. The 
theoretical equation for thermal flux was derived for an infinite medium; 
whereas, the experimental thermal flux was actually determined from 
a finite medium. The experimental phase was effected by neutron leakage, 
i1.e@., neutrons were lost at the extremities of the VISIFLUX howitzer. 

The infinite medium has no neutron leakage in that all are absorbed 
within the medium. The investigation into relaxation length established 
several theoretical relationships such as a small positive slope for 
relaxation length with temperature dependence and also that the magnitude 


was on the order of 10. With reference to Fig. 30, the relaxation 
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lengths of both bare and cadmium-covered Indium have been plotted as 
a function of temperature. The slope of the experimental values in 
contrast to the theory was found to be slightly negative. Recall 


that the flux avay from the source can be approximated by the expression 


-r/d 
Par) = Se 3 such that a smaller value of b implies a faster flux 
r 





attenuation with distance or a larger leakage loss. Since raising the 
temperature of the inedium has the effect of increasing the neutron 
kinetic energy it 1S expected that the leakage loss will increase with 
temperature. The decreasing trend of the experimental relaxation 


lengths with temperature confirms this observation. 


B. THEORETICAL VERSUS CRYSTAL DATA 

The theoretical and crystal thermal fluxes were plotted as per 
Fig. 41 for comparison with excellent agreement being obtained in 
the shape of the therinal flux curve at the temperature selected 
(T = 102° F). The general trend of the thermal flux as determined by 
the crystal also agrees with the theoretical analysis, i.e., an increase 
in flux magnitude was obtained for an increase in temperature. Crystal 
thermal flux magnitudes with distance and temperature dependence are 
tabulated in Table 8 as determined from the computer program listed 
in Appendix C. 

A minor drawback in using the crystal was the nonavailability of 
a standard for calibration of the crystal. An overall system efficiency 
of 25 per cent was assumed which appears to be realistic in view of the 
correlation of two sets of data. 

The advantages of crystal detectors for flux measurement are its 
high efficiency and its rapid acquisition of data. The high efficiency 
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of a crystal in a neutron detector permits the usage of a detector 


of a very small size which enhances the accessibility of the detector, 
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TABLE 10 


THEORETICAL THERMAL FLUX AS A FUNCTION 


OF DISTANCE, TEMPERATURE AND NON-INTEGER m 
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4694505 
3686139 
2905993 
2300053 
1827567 
1457689 
1167001 
937659 
756021 
611622 
496402 
404134 
329988 
270202 
221839 
182592 
150649 
124575 
103234 
85721 
71314 
59434 
49616 
41486 
34740 
29131 


m= .28] 


102 


4714064 
3720631 
2947335 
2343242 
1869629 
1496975 
1202706 
969505 
784048 
636046 
517535 
42232) 
345577 
283526 
232202 
192270 
158882 
131576 
109185 
90779 
75614 
63091 
52727 
44134 
36995 
31053 
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Temperature (°F) 
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1908850 
1535051 
1238321 
1062005 
813197 
661866 
540197 
442076 
362708 
298324 
245949 
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168284 
139635 
116089 
96691 
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47311 
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33386 
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4707196 
3746639 
2991607 
2396280 
1925407 
1551805 
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826933 
674214 
S572 
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305712 
232300 
208783 
1730388 
1437381 
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4681372 
3743658 
3002633 
2415365 
1948583 
157650) 
1279035 
1040549 
848776 
694150 
569121 
467747 
385333 
318157 
263262 
218291 
181363 
150969 
125896 
105170 
83000 
73749 
61898 
52025 
43785 
36398 
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TABLE 10 (CONTINUED) 


THEORETICAL THERMAL FLUX AS A FUNCTION 


OF DISTANCE, TEMPERATURE AND NON-INTEGER m 
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2446) 
20564 
109 
14585 


4617565 
3625726 
2858366 
2262357 
1797614 
1433799 
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488266 
397511 
324580 
265774 
218203 
179600 
148180 
122533 
101542 
84316 
70145 
58460 


102 


26099 
Zea 


18503 


15605 


in 


4597102 
5625318 
2837/4208 
eC oOe 
1823241 
1459833 
1172865 
945450 
764594 
620265 
004694 
411842 
337003 
27649] 
227416 
187499 
154940 
Zool 
10647 
88527 
73738 
G1525 


88 


Temperature (°F) 
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TABLE 10 (CONTINUED) 


THEORETICAL THERMAL FLUX AS A FUNCTION 


OF DISTANCE, TEMPERATURE AND NON-INTEGER m 
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48803 
40806 
34170 
28654 
24060 
20227 
17025 
14346 
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43039 
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13044 
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Temperature (°F) 
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VI. CONCLUSIONS 


The crystal and light pipe combination exhibit excellent potential 
for neutron flux determination as a function of position and temperature. 
The major advantages of the phosphor configuration used in this investi- 
gation were speed and flexibility. Data can be obtained rapidly and 
evaluated in a straigntforward manner. In contrast, the Indium foil 
method was found to be extremely time consuming because of irradiation 
time requirements, foil counting time and extensive data analysis and 
calculations required prior to nese the data. 

The Indium foil technique will undountedly provide the most reliable 
data providing all possible sources of error are minimized; however, the 
crystal combination remains consicerably simpler and much faster for 
data analysis. The crystal is adaptable to many different usages because 
of its size. In addition, the crystal would be suitable as a detector 
for isotopes with short half-lives. 

The experimental phases of this investigation verified the trend of 
increasing flux with increasing temperature as predicted by theory. 

Since the cross section generally decreases with temperature, it is 
expected cheve he neutron reaction rates are almost insensitive to 
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temperature, 
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SATURATED SAT ACTIVITY SAT ACTIVITY 
ACTIVITY TIMES R° TIMES R’ 
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COMPUTER DATA TABULATION FOR T = 102 F 
CADMIUM COVERED FOILS 
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COMPUTER PROGRAM 4 
COMPUTER DATA TABULATION FOR T = 129° F 
CADMIUM COVERED FOILS 


DISTANCE SATURATED SAT ACTIVITY SAT ACTIVITY SIGMA 
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CAD IU COVERED FOILS 


SAT ACTIVITY 
TIVES Rt 
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